Quantum Transport in Graphene Nanoribbons with Realistic Edges 
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Due to their unique electrical properties, graphene nanoribbons (GNRs) show great promise as 
the building blocks of novel electronic devices. However, these properties are strongly dependent 
on the geometry of the edges of the graphene devices. Thus far only zigzag and armchair edges 
have been extensively studied. However, several other self passivating edge reconstructions are 
possible, and were experimentally observed. Here we utilize the Nonequilibrium Green's Function 
(NEGF) technique in conjunction with tight binding methods to model quantum transport through 
armchair, zigzag, and several other self-passivated edge reconstructions. In addition we consider the 
experimentally relevant cases of mixed edges, where random combinations of possible terminations 
exist on a given GNR boundary. We find that transport through GNR's with self- passivating edge 
reconstructions is governed by the sublattice structure of the edges, in a manner similar to their 
parent zigzag or armchair configurations. Furthermore, we find that the reconstructed armchair 
GNR's have a larger band gap energy than pristine armchair edges and are more robust against 
edge disorder. These results offer novel insights into the transport in GNRs with realistic edges and 
are thus of paramount importance in the development of GNR based devices. 



I. INTRODUCTION 

The study of quantum transport in graphene nanorib- 
bons (GNR's) is of particular interest for the develop- 
ment of novel nanoelectronic devices since, unlike in 
conventional materials, the transport characteristics of 
nanostructured graphene based devices are affected by 
the geometry (chirality) of the GNR edges. In the ideal 
case the edges of GNR's can be arranged in two dis- 
tinct geometries: zigzag and armchair, as shown in Fig. 
1A. Pristine armchair edges, i.e., armchair edges with no 
atomic roughness, are semiconducting [U[2], while pris- 
tine zigzag edges are metallic [TJ [2], as shown in Figs. 
IB and 1C, in which the current density through zigzag 
and armchair terminated GNR's is plotted, as calculated 
using the methodology described below. The homoge- 
nous current density in the zigzag GNR is a signature of 
the metallic density of states at the edges. Conversely, 
the current density of the armchair terminated GNR is 
inhomogeneous and shows several distinct paths of elec- 
tronic carriers, known as quantum billiards 011], which 
result from the reflection of carriers off of the semicon- 
ducting edges. The semiconducting nature of the arm- 
chair edge originates from a quantum confinement effect, 
whereas the metallic nature of the zigzag edge stems from 
a nondispersive state localized at the periphery of the 
GNR [5 , 6j. The localized states stems from the fact that 
the bottom of the conduction band and the top of the va- 
lence band are always degenerate at the k = ir point in 
the Brilluoin zone. The origins of this degeneracy are 
discussed in detail in Ref. |2J. 

The contrasting behavior of chiral zigzag and arm- 
chair edges allows the possibility of utilizing armchair 
and zigzag GNR's for different electronic and photonic 
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Figure 1: A A schematic of a graphene lattice highlighting 
the ideal armchair and zigzag edges. The edge atoms are col- 
ored red and blue to reflect the valleys to which they belong. 
Contour plots of the 2D current density through perfect B 
zigzag and C armchair nanoribbons, for the device geometry 
shown in panel D. The reddish (bluish) colors correspond to 
high (low) current densities. Both ribbons are 10 nm x 30 
nm, backgate Vt> g bias is 0.01 eV, and an infinitesimal source 
drain bias. 



applications. For example, armchair GNR's are ideal for 
use in high speed transistor [71 18] and photodetector [9] 
devices that require a bandgap. The localized state of 
the zigzag GNR's is ferromagnetic [10] and is thus capa- 
ble of rendering the GNR half-metallic under the appli- 
cation of a lateral electric field [TT] . Consequently, zigzag 
GNR's are ideal for spintronic applications. Furthermore 
combining zigzag and armchair edges in a single device 
results in novel behaviors such as current rectification in 
Z-shaped GNR's [H[13] and spin filtering [14]. 

Although extensive theoretical models have been put 
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forth to study transport through pristine armchair and 
zigzag GNR's, [15J little is known about transport 
through GNR's with realistic edges, i.e., edges that are 
not atomically pristine. Since all experimentally re- 
alized GNR's have thus far possessed some degree of 
edge roughness [161 02] > an understanding of transport 
through these devices is vital for the design of novel elec- 
tronic devices. Furthermore, while ideal armchair and 
zigzag edges are the most frequently studied edge types, 
other self-passivating edge reconstructions are possible 
[18J and have been experimentally observed [19 , 20J. The 
influence of partially and fully reconstructed edges on 
transport in GNR's has not yet been analyzed in de- 
tail. In this work we utilize nonequilibrium Green's func- 
tions (NEGF) methods to numerically simulate quantum 
transport through GNR's with both atomically smooth 
as well as rough edges. We furthermore simulate trans- 
port in self-passivated GNR's with various degrees of 
edge roughness. We find that transport through GNR's 
with self-passivating edge reconstructions is governed by 
the sublattice structure of the edges, in a manner similar 
to their parent zigzag or armchair configurations. Fur- 
thermore, we find that the reconstructed armchair GNR's 
have a larger band gap energy than pristine armchair 
edges and are more robust against edge disorder. Taken 
together our results elucidate the effects of edge recon- 
structions on the transport properties of GNR's and are 
thus of paramount importance for the design of novel 
electronic and photonic devices. 



where Gf- (E) is the retarded nonequilibrium Green's 
function given by 

G<(E) = i J dte~* Et (^(t)). 

The systems considered in this study consist of narrow 
graphene nanoribbons, typically ~10 x 30 nm connected 
to source-drain electrodes via Ohmic contacts. The Fermi 
energy Ep may be adjusted via the application of a back- 
gate bias [3] as shown in Fig. Id. The above formalism 
is sufficient to describe transport through ideal armchair 
and zigzag GNR's. However, in order to properly account 
for the altered bond lengths of the self-passivated edge re- 
constructions [18J, the environmentally dependent tight 
binding model was utilized |22| . in which the change in 
the bond lengths is accounted for by modifying the values 
of the hopping integrals ji according to 

7, (R jk ) = cni?-r ex p (-«3i?-r) [i - s jk ] , (3) 

where are the scaling parameters and Sjk are the 
screening functions calculated from first principles cal- 
culations in Ref. [22J. Rjk are the bond lengths of the 
edge reconstructions, which are calculated in Ref. [18J. 
Following the results of Ref. [22] , the screening function 
is of the form 



A. Model 

Numerical simulations have been carried out with the 
aid of nonequilibrium Green's functions (NEGF) in con- 
junction with the tight binding formalism, where only 
the first order (nearest neighbor) hopping is considered. 
The simulations were implemented using the KNIT al- 
gorithm [21] on a 46 node parallel cluster. The system 
Hamiltonian is of the form 



(i) 



where 7^ is the nearest neighbor hopping integral, €{ is 
the on-site energy, and c\ / Cj are the creation / annihi- 
lation operators for the i-th and j-th lattice sites respec- 
tively. With the aid of the above Hamiltonian the local 
current Iij and density of states pi can be computed using 
the standard NEGF formalism according to 



J dE[ 7ij G< i (E)- lij G< j (E)] and 



1 

2^ 
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(2) 



>kj 



exp (gjfc) - exp(-gj fc ) 
ex P (€jk) + exp(-£ jk ) 



with 



ijk = ft^exp 



Rji + Rki 
Rji 



(4) 



where are the screening, and Rji and Rki are the next- 
nearest neighbor bond lengths. The bond lengths and the 
values of the nearest neighbor hopping integrals are listed 
in Fig. 2f. 

Before proceeding to the discussion of our results, it 
should be noted that spin effects were previously pre- 
dicted to play an important role [231 24J in transport 
of pristine zigzag GNRs whose edge states result in flat 
bands at the Fermi energy. Spin effects may be in- 
cluded by adding a Hubbard term H' = J/^n^n^, to 

i 

the Hamiltonian where U is a parameter that defines 
the strength of the on-site Coulomb interaction. Pre- 
vious calculations with zigzag GNRs have shown that 
i^'leads to the development of ferromagnetic spin orien- 
tations along each edge and antiferromagnetic correla- 
tions between the two edges. The spin polarization of 
the edge states breaks the degeneracy of the flat energy 
bands at the Fermi energy with a gap opening up at 
Ep. This band gap leads to a vanishing of the density of 
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Figure 2: Illustrations of the edge reconstructions and their densities of states relative to the Fermi energy for the A armchair, 
B AC677, C AC56 D zigzag, and E ZZ57 nanoribbons. The table in F shows the values of the bond lengths (adopted from 
Ref. 22) and change in the hopping energy ji/jo, where the values of 7$ (i = 1, 2 • • • 10) are indicated in the diagrams in panels 
a-e. In these simulations we used 70 = 2.8 eV. 



states at Ep while higher energy states are unaffected by 
H ' [23j [25] [26] . There is some experimental evidence for 
the energy splitting of the edge states [27- 29J although 
there is a debate as to whether such magnetism could 
exist in real world conditions [26J. Among the many ef- 
fects that would quench this edge state magnetism is edge 
reconstruction. Edge reconstructed zigzag GNRs are ex- 
pected to show no magnetic effects since the edge state 
energy bands at Ep are no longer fiat [35j [36]. Here we 
do not consider nonlinear spin effects in our results. 



II. RESULTS AND DISCUSSION 

Figure 2 shows the atomic arrangements of the five 
types of edges considered in this work: the armchair 



and zigzag edges as well as the three types of self- 
passivating reconstructions, and their corresponding den- 
sities of states. Although pristine armchair and zigzag 
edges are expected for the honeycomb lattice geometry of 
graphene, thermal [30] and optical [17] edge reconstruc- 
tion has been observed to give rise to the self passivating 
variants of the pristine edges. The most common mecha- 
nism of bond rearrangement in sp2 crystalline carbons is 
the Stone- Wales mechanism, which in the case of GNR 
edges can lead to the ZZ57 reconstruction. This mecha- 
nism has been studied extensively in relation to nanocrys- 
talline carbon systems [31-33J. Furthermore, the ener- 
getics of other self-passivating edge reconstructions were 
also investigated [18J. The names of the reconstructed 
edges originate from the parent edge type, i.e., zigzag or 
armchair, and the geometry of the closed structure. For 
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Figure 3: A Conductance of ideal armchair (solid blue), 
AC677 (dash-dotted green), and AC56 (dashed pink) lines. 
C Conductance through ideal zigzag (solid black) and ZZ57 
(dotted red) lines. The 2D current density for B AC56 and 
D ZZ57 reconstructions for the same parameters as in Fig. 1. 



example, the ZZ57 edge is a zigzag edge that is termi- 
nated by alternating pentagons and heptagons. As can 
be seen in Fig. 2, the self passivating edge reconstruc- 
tions together with the unaltered edges form two distinct 
families of edges: the zigzag and the armchair family. Al- 
though the reconstructed edges bear little resemblance 
to their parent edge types, the zigzag reconstruction is 
similar to the pristine zigzag edge in that it possesses 
a metallic state at E = 0, as can be seen in Figs. 2D 
and 2E. Interestingly, the magnitude of the E = state 
is found to be larger for the ZZ57 reconstruction than 
for the zigzag GNR of the same width. Consequently, 
ZZ57 edges might be of greater utility for devises that 
exploit the edge localized states, such as devices based on 
graphene antidot superlattices |6J. Phenomenologically, 
the increase in the magnitude of the E = state can 
be attributed to Anderson localization that arises from 
the additional break in the periodicity of the hexagonal 
lattice at the ZZ57 edge, as can be seen in Figs. 2D 
and 2E. Similarly both of the armchair reconstructions 
are semiconducting and possess a finite energy gap as 
discussed in detail below. The similarities between the 
parent edge types and their reconstructed counterparts 
can be readily understood by considering the sublattice 
structure of the terminating atoms at the edges. Both 
zigzag and ZZ57 edges are terminated by atoms that be- 
long to the same sublattice, whereas armchair as well 
as AC677 and AC56 edges contain atoms that belong to 
both sublattices. The terminating atoms are atoms that 
are the outermost atoms of the GNR. Consequently, the 
armchair edges support intervalley scattering, whereas 
the zigzag edges do not. Thus the sublattice structure 
of the terminating atoms in the self passivating edge re- 
constructions determines the properties of the resultant 
edges in a manner identical to the pristine armchair and 
zigzag edges. 

Having discussed the principle features of the zigzag 
and armchair families we now turn to the discussion of 



transport through GNR's with self passivating edge re- 
constructions. We first restrict the discussion to fully re- 
constructed edges. Figures 3A and 3C show the conduc- 
tance of pristine and reconstructed armchair and zigzag 
terminated GNR's respectively. Figures 3B and 3D show 
the 2-D current density of AC56 and ZZ57 edges. Fur- 
thermore the 2-D current density for the AC677 recon- 
struction was found to be similar to that of the AC56 
edge (data not shown). The prominent features in the 
conductance through GNRs with reconstructed edges are 
the dips in the conductance at certain values of the Fermi 
energy. These dips arise from the fact that the size 
of the unit cell is enlarged at the reconstructed edges. 
Consequently, fewer states are available per unit area 
as compared to the unreconstructed edges resulting in 
pronounced dips in the conductance. It should be fur- 
ther noted that these dips appear energetically far away 
from the Dirac point. These features can potentially be 
probed in either ARPES or optical conductivity experi- 
ments. However, since no dips are observed in the vicin- 
ity of the Dirac point, the low bias conductance of recon- 
structed GNR's should be similar to the conductance of 
pristine GNR's. Therefore, our results demonstrate that 
while transport in narrow GNR's is highly sensitive to the 
type of edge terminations, namely, zigzag or armchair, it 
is robust with regard to self-passivating edge reconstruc- 
tions within the same edge family. Consequently, the 
relevant physical mechanisms that dominate transport 
in armchair and zigzag terminated GNR's can be still 
observed in GNR's with realistic edge reconstructions. 

The 2D current density of AC56 and ZZ57 edges fur- 
ther highlights the similarity between the reconstructed 
and pristine GNR edges. The AC56 terminated GNR's 
show an oscillatory current density characteristic of quan- 
tum billiards in pristine armchair GNRs, whereas the cur- 
rent density in the ZZ57 terminated GNR is homogenous. 
The most striking difference in the 2D current density be- 
tween pristine and reconstructed edges is the presence of 
edge-localized charge density. Phenomenologically this 
effect can be understood in terms of Anderson localiza- 
tion since the reconstructed edges break the periodicity 
of the underlying lattice thereby localizing the electronic 
wavefunctions. As a result the quantum billiard patterns 
are smeared by the edge- localized modes. 

Having analyzed the effects of complete reconstruction 
on transport through GNR's we now turn our atten- 
tion to the case of partial reconstructions. To study the 
effects of partial reconstructions on transport through 
GNR's we first calculate the conductance of a pristine 
armchair and zigzag GNR and then introduce various 
percentages of AC677 and ZZ57 edges respectively. To 
quantify the effect of reconstruction, we plot the conduc- 
tance value of two dips (E — 2.1 eV for the armchair 
reconstruction and E = 1.1 eV for the zigzag reconstruc- 
tion) as shown in Fig. 4. The dips are identified by 
black arrows in Figs. 3Aand 3C. As can be seen in Fig. 
4 the onset of the dips occurs at small percentages of 
edge reconstructions (~5%). Increasing the amount of 
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% edge reconstruction 

Figure 4: Conductance of an armchair (red triangles) and a 
zigzag (green circles) nanoribbons with varying percentages 
of edge reconstruction with AC677 and ZZ57 substitutions 
respectively. The conductance values were extracted for Ef = 
2.1 and Ef = 1.1 eV for the armchair and zigzag nanoribbons 
respectively, as indicated by the arrows in Fig. 3a and 3c, for 
the same size devices as indicated in Fig. 1. 



edge reconstruction increases the amplitude of the dips, 
which become most pronounced for fully reconstructed 
GNR's. It should be noted that the armchair recon- 
structions have a more dramatic effect on the conduc- 
tivity than the zigzag reconstructions. The magnitude of 
the dip in the conductance of an armchair GNR recon- 
structed by AC677 segments decreases exponentially ac- 
cording to G (%£ A C667) = M exp (-%XAcm7/A 2 ) + A 3 , 
where %#AC667 * s ^ ne percentage of AC677 segments and 
A\ = 8.1, A2 = 44.1, and As = 13.0 are empirical fit 
parameters in units of e 2 /heV. The magnitude of the 
dip in the conductance of a zigzag GNR decreases lin- 
early as a function of ZZ57 reconstruction according to 
G(%x Z z57) = Bi (%£ Z Z57) + B2, where %x Z z57 is the 
percentage of ZZ57 reconstruction and B\ = —0.17 and 
B2 = 34.2 are empirical fit parameters in units of e 2 /heV . 

Finally, we turn our attention to the bandgap of recon- 
structed armchair GNR's. In light of the above discussion 
it is sufficient to only consider fully reconstructed GNR's 
since only the absolute magnitude of the bandgap is ef- 
fected by partially reconstructed GNR's. The gap energy, 
Eg, is extracted from the conductance of the GNR. Fig- 
ures 5a and 5b show the magnitude of Eq as a function 
of ribbon width for AC677 and AC56 ribbons respec- 
tively. As can be seen in Fig. 4, the bandgap of very 
narrow GNR's less than 4 nm surpasses 1 eV; however, 
in both cases it decreases exponentially as a function of 
ribbon width. For the AC677 terminated ribbon the gap 
energy as a function of ribbon width E^ 0667 " 1 (n) de- 
creases according to E^ C667 \n) = C\ exp (— n/C2)-\-Cs 1 
where n is the number of atoms that span the width TV 
of the GNR, and d = 2.1, C 2 = 17.5, and C 3 = 0.2 
are empirical fit parameters in units of e 2 /heV. Simi- 
larly, E^ C56) (n) for the AC56 ribbons varies according; 

to E^ C56 \n) = D x exp (-n/D 2 ) + D 3 , where D 1 = 1.7, 
D2 = 17.8, and D3 = 0.2 are empirical fit parameters in 
units of e 2 /heV '. A striking difference between the AC677 
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Figure 5: Bandgap as a function of ribbon width for A AC677 
and B AC 56 nanoribbons, for the same device parameters as 
indicated in Fig. 1. n is the number of Carbon atoms that 
span the width N of the GNR. The different symbols in panel 
'a' represent the multiplicity of the number of atoms that 
span the width of the GNR. The circles, crosses, and squares 
indicate GNR's of widths that are multiples of N = 3n, N = 
3n + 1, and N = 3n + 2 respectively. 



and AC56 GNR's is that the AC677 possess an additional 
family dependence of the gap energy on the multiplic- 
ity of the number of atoms that span the width GNR. 
Namely, AC677 terminated GNR's with width's that are 
multiples of 3n and 3n + 1 possess an electronic bandgap, 
whereas AC677 terminated GNR's with widths that are 
multiples of 3n + 2 do not have a bandgap. The same 
family dependence on the width multiplicity was previ- 
ously calculated for pristine armchair terminated GNR's 
[34H36] . This behavior, however, vanishes in AC56 ter- 
minated GNR's. The lack of the multiplicity dependence 
of Eq in AC56 terminated GNR's can be understood 
from the fact that the AC56 geometry imposes a greater 
strain on the underlying lattice, thereby distorting more 
of the features of pristine armchair GNR's. Finally, it 
is interesting to note that Eq in both AC677 and AC56 
terminated ribbons is twice as large as Eq in pristine 
armchair nanoribbons of the same width (see, for exam- 
ple, the numerical calculations of the armchair bandgap 
in Ref. [34]). Consequently, these variants of the arm- 
chair GNR's might become more useful candidates for 
GNR based devices, if a reliable method for fabricating 
these edges is found. 
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III. SUMMARY 

In summary, we simulated quantum transport through 
graphene nanoribbons with realistic edge terminations 
by utilizing nonequilibrium Green's functions. We found 
that transport in the self-passivated edge reconstructions 
is qualitatively similar to transport in the parent edge 
type. This similarity is attributed to the sublattice struc- 
ture of the terminating atoms at the edges. We also found 
that all of the self passivating reconstructions localize 
electron wavefunctions at the edge. This localization may 
be understood in terms of Anderson localization since the 
reconstructed edges further break the periodicity of the 



underlying lattice. We further analyzed the width de- 
pendence of the bandgap in the reconstructed armchair 
GNR's, and found that the gap energy is twice as large as 
the gap energy of pristine armchair GNR's of the same 
width. Finally, we have shown that the band gap in 
AC56 terminated GNR's is more robust against the mul- 
tiplicity, which is of advantage for fabrication which in- 
evitably introduces width fluctuations. Taken together, 
these results elucidate the nature of quantum transport 
in graphene based devices with realistic edges and are 
thus important for the development of future graphene 
based nanoelectronic devices. 
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